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Abstract 

We derive the bosonization rules for free fermions on a half-line with phys- 
ically sensible boundary conditions for Luttinger fermions. We use path- 
integral methods to calculate the bosonized fermionic currents on the half- 
line and derive their commutation relations for a system with a boundary. 
We compute the fermion determinant of the fermionic fluctuations for a sys- 
tem with a boundary using Forman's approach. We find that the degrees 
of freedom induced at the boundary do not to modify the commutation re- 
lations of the bulk. We give an explicit derivation of the bosonization rules 
for the fermion operators for a system with boundaries. We derive a set of 
bosonization rules for the Fermi operators which include the explicit effect 
of the boundaries and of boundary degrees of freedom. As a byproduct, we 
calculate the one-particle Green's function and determine the effects of the 
boundaries on its analytic structure. 

PACS numbers: 11.10.Kk,72.15.Qm, 73.20.Dx 
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I. INTRODUCTION 



The Fermi-Bose equivalence, commonly known as bosonization, is one of the most power- 
ful approaches to extract the non-perturbative behavior of 1 + 1-dimensional Field theories 
and strongly interacting Condensed Matter systems. This approach is an outgrowth of the 
work on the solution of the Luttinger-Thirring model by a long list of authors The 
Fermi-Bose equivalence is deeply rooted in the properties of the algebras of the currents 
and densities of the physical observables. An alternative picture of the same physics is also 
found in the functional integral approach |7],[|]. In its standard form, bosonization maps a 
system of interacting fermions into a system of bosons. The system of fermions (as well as 
the equivalent bosonic theory) is usually defined on manifolds without boundaries. This is 
appropriate for the study of the thermodynamic limit. In the case of the simple Luttinger- 
Thirring model, which has an abelian U(l) symmetry, the bosonic theory is essentially free 
and the solution of the interacting fermion theory is thus found. For more complicated situa- 
tions, in which the symmetry is non-abelian, the bosonized theory is a non-trivial Conformal 
Field Theory [|-[12[ . 



There are a number of situation of physical interest in which systems with boundaries 
are important. Of particular interest are the so-called quantum wires. Also, there are 
a number of situations in which a localized degree of freedom interacts with an extended 
system of fermions. Examples are one- dimensional systems ( as for instance, quantum 
wires) interacting with magnetic and non-magnetic impurities. The interplay of Luttinger- 
like behavior and the Kondo effect is a subject of considerable interest fl3| . |l4}l . Similar type 



of physics is also found in the context of the Callan-Rubakov effect [15]. A number of authors 



have used bosonization methods to the study of semi-infinite systems JTj] and systems with 



quantum impurities [16] such as the Kondo effect [17] and spin chains with impurities [18 



However, in all cases, the bosonization rules of the bulk system were used and little attention 
was paid to the role that the boundary conditions of the fermions may play in the bosonized 
theory. 
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In this paper we discuss the effects of the boundary conditions of the fermions and in 
which way they affect the physics of the bosonized theory Since bosonization is seen most 
clearly in the case of free fermions, we consider here just the simplest case, namely non- 
interacting, abelian, theory of fermions on a half-line ("~ + 1-dimensions" ) . In a separate 



publication [19] we will discuss the case of interactions and the Kondo effect, which involves 
systems with non-abelian symmetries. We consider systems of fermions with two types of 
boundary conditions: (a) R = — L at the boundary and (b) R(0,X2) = — e 2ld ^L(0, x 2 ), 
again at the boundary. This choice of boundary conditions is motivated by the underlying 
microscopic physics of the systems of interest. In the case of quantum wires, the microscopic 
system consists of non-relativistic fermions at finite density. At a physical edge (i. e. the 
boundary) the fermion amplitude ^(0) must vanish. Upon a decomposition into Left and 
Right movers, we find that a sharp edge implies the choice (a), i. e. R{0) = —L(0). If a 
localized quantum impurity, parametrized by a degree of freedom 9, is placed at the edge, 
the boundary condition changes to case (b),z. e. R{0,x 2 ) = -e 2i9 ^L(0,x 2 ). Physically, 
this boundary condition is caused by a coupling between the fermion charge density at the 
edge and the localized quantum impurity. We also introduce physically sensible models of 
the quantum dynamics of the impurity. 

In section |T|, we study a system of free electrons confined to the half-line and derive an 
expression for its partition function. In order to write this partition function only in terms of 
bosonic operators, in Section |TXT| we compute the determinant of the 1 + ^-dimensional Dirac 



operator coupled to a source which is defined as a delta function at the boundary. We show 
that this problem is equivalent to the computation of the determinant of the Dirac operator 
plus a suitable boundary condition . This last determinant is solved by using Forman's 
approach [pH| . 

In section |TV], we derive an expression for the partition function only in terms of bosonic 
degrees of freedom. In other terms, we present an alternative form of abelian bosonization in 
terms of functional integrals. Although this is a subject that has been discussed extensiveley 
in the literature, we present here a derivation of this classic result in a form that is suitable for 



treating systems with boundaries. We use these results to calculate the bosonized currents, 
the current correlation functions and their anti-commutation relations paying close attention 
to the role of the boundary conditions of the fermions. 

In section [V], we derive the one particle fermionic Green's function from the bosonized 
theory. We discuss two ) when the boundary condition is R = —L , and (b) when 

there is an extra degree of freedom at the boundary, i. e., R(0, X2) = —e 2l6 ^ X2 'L{Q 1 X2)- We 
focus on the way the one particle Green's function is affected by the presence of the boundary 



and of a dynamical degree of freedom at the boundary. In section [Vj we discuss a model of 
free fermions coupled to a dynamical boundary degree of freedom, and calculate the fermion 
one-particle Green's function for this case of interest. We consider two asymptotic regimes: 
(a) one in which the boundary degree of freedom is strongly pinned, and (b) the case in 
which it fluctuates wildly. In the Conclusions we give a brief summary of our results. In the 
appendices we give details of the computation of the boundary contributions to the fermion 
Jacobian, of the use of Forman's method and to Euclidean-Minkowski correspondences are 
given in the Appendices. 



II. FREE FERMIONS ON A HALF-LINE 

In this section we study a system of free massless Dirac fermions confined to the half 
line (i. e. , 1 + ~ dimensions) satisfying dynamical boundary conditions. Our goal is to 
compute the partition function in order to be able to derive the bosonizations rules for 
this system. Throughout, we work in Euclidean space. In practice, this means that the 
space-time manifold is a semi-infinite long cylinder, with its axis along the space direction 
x\ and periodic (antiperiodic for fermi fields) boundary conditions along the imaginary time 
direction x^ of perimeter T, with T — > 00. 

The (Euclidean) Lagrangian for the system (coupled to sources) is 

£, _ r r- _4 / ,./ / , (2.1) 
where is a source (a background gauge field) which couples to the fermion current = 

4 



ip "f^ip- In a Fermi theory in Euclidean space, not just the time coordinate has to be continued 
analytically to the imaginary axis. The structure of the 7-matrices and the requirement that 
the Euclidean Dirac operator be hermitean forces an analytic continuation of the chiral angle. 
Hence, the chiral transformations generated by 75 become complexified in Euclidean space. 
However, we must keep in mind that these two analytic continuations are logically distinct. 
Moreover, these two analytic continuations need to be done in order to have a physically 
sensible interpretation of the results and to recover results back in Minkowski space. Hence, 
the fermion boundary conditions themselves also have to be continued analytically. We thus 
demand that the fermions satisfy the boundary condition R(0, X2) = —e 29 ( X2 ^L(0,X2)- 
The dynamics of the fermions is described by the functional integral 

Z[A] = J V$T>ipexp (- J d 2 xC F ^j . (2.2) 

In what follows we will use the path integral bosonization method based on Seeley's 
expansions for complex powers of elliptic operators pifl . 

The first step in the bosonization process is to decouple the fermions from the sources 
by means of a combination of (suitably chosen) gauge and chiral smooth, single-valued 
transformations of the form 

$(x) = x(x) e - ir ' (x)+ ^ ix) , (2.3) 

provided that the vector potentials A^(x) can be written in the form 

A^x) = dfl(x) - e Mi A0(z). (2.4) 

This requirement is satisfied by all topologically trivial configurations of the fields A^. The 
usefulness of the transformations of Eq. ( |2.3j ) and Eq. ( |2.4j ) is that they completely decouple 
the fermions from the vector potentials. However, this transformation changes the fermionic 
measure in a non-trivial way. Namely, the change in the fermionic measure under this 
transformation is: 
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Vijj Vip = J F V\ T>x- 



(2.5) 



The jacobian, Jp , can be written as 



Jf 



Det(^ + 4) = _ j Q i w r {t)dt 
Det(^) 



(2.6) 



with 



w' — 2 / dx2dxitrK 2 [ipt'i x, x]7 5 0(a;) 



(2.7) 



where ip t = t [(l - t)A\. The kernel K 2 [ip t',x,x] is the constant term in the analytic 
expansion of the Heat Kernel < xle'^^^x >. In the case of the fermions on the full line 
(which can be considered as a close manifold if we impose vanishing boundary conditions 
at the infinity) the jacobian is well known and it has been calculated by several authors 
In the case of a manifold with boundary we have to evaluate this Heat Kernel, following 
Atiyah-Patodi-Singer (APS) ||22|| , separating the contributions of the boundary from the 
ones of the bulk. Hence, the decoupling of the gauge field has to be done in two steps, first 
we decouple the bulk and then the boundary. For this pourpose it is convenient to write the 
field as A^ = + s M where s M is defined as the value of A^ at the boundary, and as 
its value in the bulk. That is, s M = A^5{x{) and is zero at the origen. 

We define a chiral and a gauge transformation as in Eq. (|2.3|) , such that it only decouples 
the field B^ from the fermions. Therefore at this point we are going to assume that the field 
B^ is the one that can be written in terms of (ft and r\ as in Eq. ( |2.4| ) . The partition function 
becomes 



where the integral in the action is restricted to the half-plane Q = {x\ X\ > 0}. The Jacobian 
reads, 




(2.8) 



Jf 




(2.9) 
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The APS procedure to compute the Heat Kernel consists in dividing the domain (Q) in two 
intervals, (0, e) and (e, oo). In the interval (e, oo) we can compute w'{t) in a usual way |7J 
and obtain the well-known result for closed compact manifolds 

where f2(e) = {x\ x% > e}. The contribution of the boundary to the jacobian is given by the 
regular term of the following expresion, (in the limit s — > 0) 

5J F = exp(^J dt J dx 2 J dx 1 < a;|e~ s ^ )t ^ )f \x > cf)(x 1 , x 2 )) . ( 2 -H) 

We show in Appendix A that this contribution is given by 




(2.12) 



Therefore, if we choose <j)(0, x 2 ) = the contribution of the boundary term to the jacobian 
is one. It is important to remark that this cancellation happens only because we made a 
chiral transformation that is trivial at the boundary. The price for this choice is that we 
still have the fermions coupled to the field s M with support only on the boundary. We will 
show how to deal with this problem in the next section. 

Then after the chiral transformation is performed the partition function is 

Z[B^ s M ] = J VxD X exp (- ^ d 2 x[ X ^+ fa]) ™P ("^ J^W? + \ J Q ^ s A<j>) 

(2.13) 

where we took the limit e — > 0. 
Since 

we can write 

d u d u (f)(x) = e a ^d a B^(x) (2.15) 
7 



(2.14) 



with the boundary condition 

0(O,x 2 ) = O (2.16) 

Therefore 

<f>(x) = J^d 2 y G (x,y)e afl d^B^y) (2.17) 

where Go(x,y) is the Green's funcion defined as 

d 2 G (x,y) = S(x-y) 

(2.18) 

G (0 1 x 2 ;y) = 0. 

Since is only defined for x% > 0, we have to use the Green's function for the half line 
(Eq. ( |2.18|) ) which is given by 

G (x, y) = g(x 2 - y 2 , xi - Vi) - g(x 2 - y 2 , x\ + yi) (2.19) 

where 



1 , (u 2 + v 2 + a 2 ) 

g(u, v) = —\n K - >- (2.20) 

Air a z 



with a acting; Kulator. 

Using these expressions to write <f>(x) in terms of , the jacobian J F can be written as 



lnJ F = -— f (d^4>) 2 + — [ e^SfJd^ 
In Jn re Jn 



= — [ dx f dy(B^(x)T flu (x,y)B u (y) + 2B^(x)T tlu (x,y)s u (y)) 

= — dx dy{B tl {x) + s fl (x))T^(x,y)(B u (x) + s v (x)) 
ZTi Jn Jn 

-77" / dx f dys fM (x)T IJil/ (x,y)s l/ (y) (2.21) 
z7r Jn Jn 

where T^x, y) is defined as 

T, u (x, y) = (d^d^5, v - dMdP) G Q (x, y) (2.22) 

We can linearize the cuadratic term in B^(x) + s Al (x) by introducing a bosonic field uj 
defined in fl such that u>(0, x 2 ) = 0. It is easy to see that Eq. fl2.21|) can be expressed as 



In J F = fC J Vu exp (—^ J^^ 2 + \ j n e iw®i> u ( s n + B v) 

6XP ( _ 27r In dX j n dys ^ T ^ x ^ Sv ^) ( 2 ' 23 ) 

Since by definition, s^{x) = s^(x 2 )S(xi), we need the value of T(x,y) at X\= y\ = 0. It is 
easy to check that r i2 | Xl ,y 1 =o = r2i| a!1 ,y 1 =o = 0. Similarly, by using the expression for the 
Green's function, Eqs. (|2.19( - |2.18|) , and for the kernels T^x, y), Eq. (|2.22j) , we also find that, 



in the limit a — > 0, rii| X:l yi= o = 0. In the same way, we find that the value of T 2 2{x, y)\ xl ,yi=o 
is given by 

r »(*.»)i--«- = -; %-»)■ + „■ + ^ - *> (2 ' 24) 

where P denotes the principal value. From now on we will assume a principal value in this 
type of expressions. Then the partition function becomes 

Z[b^s^] =JCJ VxVxVu exp (- J^d 2 x[x(0+ ff)x]) 
exp (-— f {d^ujf + - f e^d^S/t + B )) 

^(-g/ (,,( £:y )- (2 ' 25) 

Now we have solved completely the decoupling of the gauge field B^. But we still have a 
fermionic integration left to perform. This integration will be done in the next section. 



III. FREE FERMIONS IN A 5(Xi)-POTENTIAL 

In this section we calculate the determinant of the (1 + |) -dimensional Dirac operator 
coupled to a S(xi) potential. Hence we have to consider the Dirac operator, [7 M (ic? M + s^)], 
where s M = s fl (x2)S(xi). The formal way to compute the determinant of this operator is by 
computing its eigenvalues. First, it is useful to regularize the 5 function as the limit e — > 

£+- 

of the function v e (x) where v e (x) = 1. It is clear that when e, 5 — * 0, the solutions of 

6 2 

[r(id, + sjty = Xip (3.1) 
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in x\ > will be the same as the solutions of the equation 



(3.2) 



plus a suitable chosen boundary condition. To find such a boundary condition we observe 
that integrating Eq. (|3.1| ) over the interval (e — |,e + |) and assuming that the field ip is 
finite in such interval, we obtain for small 5 



- 7iOA(e + ) - (71*1(^2) + 72^2(^2)) 

Such a field can be found if we look at solutions of 



dx\ i/j(xi)v e (xi) = 0. (3.3) 



diip = (7552(^2) +is 1 (x 2 ))v e (x 1 )ip(x 1 ). 



(3.4) 



Explicitly, 



ip(x) = P X1 exp[(7 5 s 2 (x 2 ) + is 1 (x 2 )) / dyv e (y)]^(e 



(3.5) 



where P X1 is the spatial ordering operator. Since the operator in the exponent of Eq. (|3.5|) 
commutes at spatially separated points because of the shape of v e (x), we may set P x = 1, 
and taking the limit 5 — > we get 



V>(e"). 



(3.6) 



It is clear now that, in the limit e — > 0, the boundary condition we are looking for becomes 



(3.7) 



Recalling that at e the fields also satisfy Be ^^^(e ) = with B the matrix given by 



B 



(3.8) 



V 11 / 



then, in the limit e — * 0, the boundary condition for the field ip is 



BU-\x 2 )^\ Xl=0 = 



(3.9) 
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where lA(x 2 ) = e^ h ^ X2 > +S2 ^ X2 >>K This equation tells us that the source s 2 (x 2 ) can ^ e though 
as a dynamical degree of freedom at the boundary that interacts with the fermions through 
the above boundary condition. In other words, the determinant we have to compute is 
det('f IJ '(id IJ )) BU ^ 1 where 

r y^(id u ) with ib such that 

V " (3.10) 

BU-\x 2 W\ Xl=0 = 0. 

The solution for this kind of determinants was first given by Forman |2(J. His theorem 
relates the determinant of the above differential operator with the determinant of an (infinity 
dimensional) matrix, 

d 



M) 



BU- 



■ ! dsTr(^-H(s)H-\s)) 
J ds 



The matrix if is a functional of the function 9{x2) + £2(2:2) • 

In order to apply Forman's methods we introduce an auxiliar parameter r such that 



(3.11) 



U(t) 



(3.12) 



where we have defined a(x 2 ) = 0(x 2 ) + s^ixz). This method relies on the knowledge of the 
space of solutions of the homogeneous equation (ij^d^ip = 0. We must consider a com- 
plete, but otherwise arbitrary, system of functions in the kernel of this differential operator, 
which do not need to satisfy any particular boundary condition. In order to simplify the 
calculations, we choose the basis that expands the kernel of (ij^d^ip = 0, satisfying 



T T 

il> n (x u -—) = -$ n {x x ,—) , 

Bij n (x 2 ,0) = Bij n (x 2 ,L) = 0. 



(3.13) 



Later on we will send L, T — > 00 . The basis is 



(1 \ 



v 1 / 



, n G Z 



(3.14) 



where w n 



(2n+l)7r 
T 



. Forman's Theorem [EDI relates the problem of evaluating the functional 



determinant of a differential operator with a parameter-dependent boundary condition (such 
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as the one defined by Eq. ( |3.10D ), to that of obtaining the determinat of an operator $ M (r), 
acting on functions defined on the boundary : 

Here the determinants on the left-hand side are defined through the ( -function regulariza- 
tion, while the right-hand side is a well-defined quantity that can be evaluated using any 
basis at the boundary The operator $ At (r) is defined as follows. We first define the Poisson 
map Pb- If the boundary conditions are defined with the operator B, then Pb is such that 
the unique solution of 

Of(x) = 0xen 

(3.16) 

Bf(x) = h(x) iGffl 

where O is an arbitrary differential operator, is given by f(x) = Psh(x). If A is another 
boundary condition, then &ab = APb- Even though Pb is difficult to evaluate, <&ab is not 
because of its main property 

Aip = $abBiP (3.17) 

for any solution of (Dip = 0. In our case $^(t) satisfies 

Bl(-\fi + T)4> = <5>^T)BU-\n)ij (3.18) 

for any ip that obeys the free Dirac equation. Let us define h™(x2) as 

h;(x 2 ) = BU-\fi)^2) (3.19) 

where Tp n {x2) = "^(O, £2). They satisfy h r p L+T (x2) = $ M (r)/i"(a;2). If we were able to write 
h r p l (x2) and /i™ +r (2 2 ) in terms of a complete set of functions in [0, L] we would have the 
(infinite) matrix ^(t) expressed in this basis, and by using conventional methods we could 
compute its determinant. The basis will be ip n (L } X2) — e tWnX2 with w n = ( 2 "+ 1 ) 7r ( n ^ z). 
It is easy to compute h™(x2) and h^ +T (x2) in such basis. Following the definition, 
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hfa 2 ) = cosh(w n L + fxa)e iWnX2 , 
K+t( x i) = cosh(w n L + (/i + r)a)e l 

In the limit L — > oo, tanh(w n L) goes to sign(u>„), where 



sign(w n ) = sign(2n + 1) 



1, n > 
-1, n < 



Then, up to a normalization constant, we can write 



(3.20) 



(3.21) 



hfa 2 ) = e 



sign(w n )/ia(x2) e iw n X2 



(3.22) 



Note that the m-th component of h™(x 2 ) in the basis under consideration is 

1 r- 
2T J-f. 

We can think of h^ m (x 2 ) as matrix elements, in the space spanned by the basis tjj n (0,x 2 ), 
of a matrix h{r). In this way, 



< i> n \h{ T )\i'm >= h 



n,m 



(3.24) 



and the operator $ M (r) can be written as 



h(jM + T)h- 1 (T)=%(T) 



(3.25) 



Then 



^ In det S„(t) = Tr^M/i + r)^ 1 ^ + r)] - j-^h' 1 {^\)\. (3.26) 



We have reduced the computation of -^ln Det&^r) to the computation of a quantity of the 
form 



Tr^hMh- 1 ^)] 

at /i — 1, where /x is an arbitrary parameter, and is a matrix of the form 



(3.27) 



hQi) 



{ [e-1 + _ [e1 ++ J 
13 



(3.28) 



(see Appendix B for notation). At this point it is important to make a conection between 
the matrix h(fi) and the so-called Toeplitz matrices. A Toeplitz matrix of order N is the 
name for a matrix of the form 

TjM = c m . n m, n = 0, N - 1 (3.29) 

where c p are arbitrary complex numbers (p = 0, ±1, ±(N — 1)). Let us denote D N as 
its determinant. Giving a set of complex numbers c p , p > 0, and under suitable conditions 



Szego 's Theorem |23| gives an expresion for 



Clearly fi N acts as a regulator. If 



lira (3.30) 



* 2tt Jo 

and 

1 r 2n 



— r dt e~ ipt C{e tt ) (3.31) 



/i = exp{— / ^C(e 1 *)} (3.32) 
2n Jo 



g p = ^- r dt e~ ipt lnC{e lt ) (3.33) 



and 



2tt Jo 

then, Szego 's Theorem states that 

jW^f) = ex P{E W-P^}- ( 3 - 34 ) 

It is obvious that [e~ a ] and [e a ] ++ are both Toeplitz blocks in the sense that they satisfy 

the definition of a Toeplitz matrix. Notice that if h(/j,) were made up by only one of these 
blocks, the determinant (6.19) would be exactly the one given by Szego' s Theorem. We are 
going to show that, even though the off-diagonal blocks are non- vanishing, and the indices 
that label the matrix entries are integers numbers, Szego' s Theorem is still valid up to an 
overall factor of 2. We can also see that Forman's Theorem provides a natural cancelation 
for the regularization factors in the Szego's Theorem. 
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In Appendix B we show that 

Tr[j^h{ii)hr x {ii)\ = tr {[«]++ - [a] — 

+ N + -(-[e Q 1_Je— ]_ + + [e a -]„_[e Q -]_ + [ e -"-] ++ [e—] ++ ) 

+ M- + ([e- Q +] ++ [e a+ ] + _ - [e- a +] + 4e- a+ ] + Je a+ L-[e a+ ]--)} (3-35) 

After anti-transform Fourier the above expression, we see that the first two terms cancell 
with each other. The third term becomes 



Tr{[a] + 4e a -]__[e- a -]_ + } = £ mi„_„|f |.._..| ( 



n-p\y \p-q[ e }q-n 
n>0p,q<0 



p-l 



a 



n+p[ c J_ p+ gF i-q-n 



EEE 

n>0p>0 q=l 

Hi Z i n>0j»0c?=l 

r e a-(z 2 ) * e -a_(z 3 ) 

[d*i]a(zi) f[dz 2 h r f[dz 3 }- 



(z 2 - Zi) J (z 3 - z x ) 2 
[d Zl }a + { Zl )e a -^\-d Zl a^z{))e- a - { ^ 

= — J dx 2 a + (x2)d 2 a^(x2). (3.36) 

Folowing the same procedure it is easy to see that the fourth term gives exactly the same 

contribution. The two terms with [a]_+ contribute with the same kind of expresion but 

with " + " interchanged with " — " and with opposite sign. Therefore 

d — 1 r 

Tr[— /i(/x)/z _1 (s/i)] = — r dx 2 [a,+{x 2 )d 2 ai-{x2) - a^(x2)d 2 a + (x 2 )} 
d[i iri J 

1 f^ W-°<5 . (3.37) 



2tt 2 J ' (x 2 - y 2 ) 2 
Eq. ( |3. 37|) coincides up to a factor of 2 with the expresion given by the Szego's Theorem for 



each one of the diagonal blocks. The regularization factors for the determinant are [at]++ 

and [a] and we see that they cancell with each other in a natural way. This expression 

coincides with the result obtained by Falomir et al. |M| , using a perturbative approach. Note 
that our result does not involve any approximation. By inserting this result in Eq. ( J3.15 ) 
and recalling that a = 9 + S2 we obtain 
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d_ ln I Vet(0) BU -i (ll+T) 
dn I Det(^) BW _ 1((U) 



^=0^=1 



2vr 2 



dx 2 dy^ 



[(9 + s 2 )(y 2 )-(6 + s 2 )(x 2 )f 
(x 2 - y 2 ) 2 



(3.38) 



Following Falomir et al. [2^] we obtain the value of the ratio of the determinants which is 
given by 



In 



Detfzc 



>BU- 



Detm B 



4tt 5 



dx 2 dy; 



[(e + s 2 )(y 2 )-(9 + s 2 )(x 2 )f 
(x 2 - y 2 ) 2 



(3.39) 



This kind of non-local expressions have been used by Anderson and Nozieres and 
Dominicis in the X-ray absorption and emission in metals [ES|and in the context of the 
Kondo problem [27]. We have then computed the last ingredient we needed in Eq ( |2.25|) 
to complete the bosonization scheme. In the next section we describe the bosonized theory 
starting from that expression for the partition function. 



IV. BOSONIZED THEORY 



Once the fermionic generating functional is known, any correlation function can be ob- 
tained by functional diferentiation. In particular, vacuum expectation values of arbitrary 
products of fermion currents 

N N 



< n j ^ x i) >=< n ^i. x i)i^{ x i) > 



(4.1) 



i=i 



i=l 



can be obtained from the fermionic partition function 



Z[A) = / V^V^exp (- / d?xt/j (i0+ A)ip 



(4.2) 



by using the regularization prescription for defining the fermion determinant, as follows 

1 SZ[A] 



<J ^ X)>= Z[A]5A,(x) 



(4.3) 



A=0 



Since, as we have shown in section II, we have to split the external field A„ into two fields, 



and s M , we define two currents, one in the bulk 



1 5Z[b^,s^ 



Z[b^ s M ] 5B ll (x) 



(4.4) 



B,s=0 
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for x\ > and the other at the boundary 



< J^0,x 2 ) >- 



(4.5) 



S,s=0 



By using Eq. ( |2.25| ) and Eq. Q3.39| ), the partition function is 



- j e^d^Sft + B tli 



cxp 



{S2(V2) - S 2 (X 2 )) 2 

, 4vr 2 J L (a; 2 - y 2 ) 2 

[(g + a 2 )(y 2 )-(g + a2)(x 2 )] 2 < 
(x 2 - y 2 ) 2 



(4.6) 



which can be written as 



e^dui^{Sf, + B u 



exp (aJr / dX2dy2e ^ X2)K ^ X2 ' y2)d ^ y2 

exp dx 2 dy 2 6(x 2 )K(x 2 ,y 2 )s 2 (y 2 ) 



(4.7) 



The kernel K(x 2 ,y 2 ) is defined as 



K(x 2 ,y 2 ) 



1 



+ -<5(x 2 - Jfe) J ■ 



(4.8) 



v 7T (x 2 -y 2 ) 2 a 

Note that this action is already linear in the external fields. Then the currents can be easly 
read from the partition function, 



JJx) = -e^duU ,for xi > 0, 

71 



(4.9) 



which is the usual bosonization formula for the fermionic current in the bulk (in Euclidean 
space) . The value of the current at the boundary is 

5u 



J^(0,x 2 ) = -t^ v d v uj\ Xl=0 + — [ dy 2 K(x 2 ,y 2 )6^{y 2 ). 

71 71 J 

Let us now compute the current-current correlation functions. They are given by 

1 5 2 \nZ[b^s,} 



(4.10) 



< J„(x)J u (y) >\ xuyi>0 



Z[bfj,, s M ] 6B^(x 2 )5B v (y 2 ) 



(4.11) 



B=0 
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< J^x)J v (y) >\ 



1 SHnZib^s,] 



yi=o 



Z[b^ s M ] 8s ll (x 2 )8s v (y 2 ) 



(4.12) 



s=0 



< Mx)J v {y) >L 1>0)yi=0 



1 8 2 hiZ[b^s^ 



(4.13) 



B,s=0 



Z[b^, s M ] 8B^(x 2 )8s v {y 2 ) 

Before computing the functional differentiation we have to integrate the bosonic fields. In- 
stead of using Eq. (|4.7| ) for the partition function we go a step backwards and recall that 



the integral over w comes from Eq. ( |2.21|) as an alternative way of writing the Jacobian Jp 
(Eq. ( |2.10| )). By using Eqs. ( |2.13| ) and ( |3.39| ) we can also express the partition function as 



Z[B^ s M ] = exp (— J (B^x) + s^{x))T^(x, y){B u (y) + s^x)) 
exp (—- J dx 2 dy 2 6(x 2 )K(x 2 ,y 2 )6(y 2/ 
exp f-J dx 2 dy 2 d(x 2 )K(x 2 ,y 2 )s 2 (y 2 ) 



Then, the current-current correlation functions are 



(4.14) 



< Jp{x)J v {y) >\ xim>0 = -T^(x,y) 



(4.15) 



< J^(x)J v (y) >| = -8^8 u , 2 K(x 2 ,y 2 ) 

,J 7T 



(4.16) 



< Jn(x)J v (y) >\ xi>0m=0 = -T^(x,y)\ Xl=0 



(4.17) 



We follow now the method suggested in reference [0] to obtain the commutation relations 
from the current- current correlation functions at equal time. 
For x, y in the bulk, we find 

1 



[Mx\J x (y)\ = — (2T 21 (x,y)) 
Ztt 



X2=V2 
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-di8(x 1 - yi) 



(4.18) 



which is the usual commutation relation (up to a factor of i which is absent in Euclidean 
space). Close to the boundary, the only non- vanishing commutator is 



[J 2 (x),J 1 (0 : y 2 )} = -(2T 21 (x,y) 

71 



) 

£2=£/2,J/l=0 

-d 1 5(x 1 ) (4.19) 
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since at equal time T 2 i(x,y)\ yi=0 = — \di5{x-i). Hence, we see that the presence of the 
boundary changes the currents but it does not change the commutation relations. 

In next section we will see how the boundary and its degree of freedom modify the one- 
particle Green's function. 



V. BOSONIZATION OF THE FERMI OPERATORS 

In this section we derive a set of bosonization rules for the fermion operators for systems 
with two types of boundary conditions: a) fixed R(0, x 2 ) = —L(0,x 2 ) and b) dynamical 
boundary conditions R(0,x 2 ) = — exp(29(x 2 ))L(0,x 2 ), where 9(x 2 ) is a boundary degree 
of freedom. Our startegy is to first calculate the one-particle Green's function in 1 + | 
dimensions with both types of boundary conditions and to use these results to determine 
the bosonization rules. We will see that we can obtain a generalization of the Mandelstam 
bosonization rules for a semi-infinite system with specific boundary conditions. In other 
words, we are going to find bosonic variables, u(x) defined in the bulk, and 9(x 2 ) defined 
at the boundary, with an action >Sb(cj, 9), such that in the framework of this bosonic theory 
there exist operators with the same expectation values that those of the fermionic operators, 
R and L in the framework of the fermionic theory. 

A. Boundary Condition R=-L 

The Fermionic one particle Green's function is 

S a f3(x,y) =< ^ a {x)^p{y) > . (5.1) 
and it satisfies the Dirac equation 
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i$ x S (x,y) = 5(x - y) 
S 1 1 (0,x 2 ;y) = -S 2 \0,x 2 ;y) 

This boundary condition implies precisely R(0,x 2 ) = —L(0,x 2 ). The solution for this dif- 
ferential equation is 

S (x,y) = S F (x,y) - >y 2 S F (x*,y) (5.3) 

where S F is the free fermionic Green's function, and x* = (— xi, x 2 ). Explicitly, 

( i i \ 



So(x,y) = 



(5.4) 



V 2 w J 

with z = (xi - yi) + i(x 2 - y 2 ) an d ™ = - Di) + ^(^2 - 2/2)- 

We will now proceed to find a bosonized expression for the fermion operator consistent 
with the boundary condition R = —L. The method that we follow is a generalization of 
the procedure used by Fradkin and Kadanoff (FK) in two-dimensional critical phenomena 
p8| , |29[| . It consists of finding a set of operators in a theory of a free bose field u> whose 
correlation function are given by Eq. (|5.4j ). Following FK, we define a bosonic action coupled 
to a background gauge field as 

S B (u, A^) = -i- jf d 2 x (d^u + A^f (5.5) 

where a; is a bosonic variable defined in Q which satisfies vanishing boundary conditions. 
Here A^ is a vector field such that 

B(z) = e^A^z) = n(-5(z - x) + 5{z - y)). (5.6) 

The gauge field A^ represents two flux tubes, each of flux tc and — ir respectively, at points 
x and y. In other terms, we have two disorder operators. Hereafter, we will denote the 
insertion at a; of a disorder operator of flux (or vorticity) ±7r by the symbol JC± n (x). For 
infinite systems, FK found that the fermion operators are represented by suitable chosen 
linear combinations of products of disorder operators JC±^(x) and order (or charge) operators 
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exp((±2u;(V)) (for x — > x'). This construction is the euclidean analog of the Mandelstam 
operators ||] of abelian bosonization. We will now investigate how does the presence of a 
boundary (as well as specific boundary conditions) alter these rules. 

We are interested in computing the expectation values of operators of the form 
e {±iu(x)±iuj{y)) - n ^ e presence of disorder operators (i.e. gauge fields of the from of Eq. ( |5.6| )). 
Therefore we need to compute bosonic partition functions of the form 

Z[A„, J} = Jvlo exp (sB + ij^ d 2 z J(z)uo(z^j (5.7) 

where the source J has the form 

J(z) = a5(z - x) + b5(z - y) (5.8) 

and a, b take the values ±1. We will show that the expectation values of these operators in 
the bosonic theory give the components of the one particle fermionic Green's function Q5.4|) . 
Since the action above is quadratic in u it can be integrated explicitly 

Z[^,J] = /Cexp {-^j n d 2 x A 2 ^ 

expf-^f d 2 xd 2 y (-d fl A fl + U)(x)G (x,y)(-d^ + iJ)(y)) (5.9) 

where the Green's function G (x,y) was defined in Eq. ( [2.18] ). Note that, as it has been 
discused before, this Green's function is consistent with the boundary condition imposed on 
uj. After some algebra (notice that the integration by parts can be safely done because of 
the boundary conditions chosen) we get for Z\A^ J] 

Z\A^ J] = K. exp jf d 2 x d 2 y A^(x)T^(x, y)A v (y) s j 

exp Q Jj 2 x d 2 y [A^x)d^G (x,y)J(y) + J{x)d^G Q {x,y)A v {y)\) 

exp f| J^d 2 x d 2 y J(x)G (x,y)J(y)^J (5.10) 

where the function T fiu (x, y) is defined as 

T Xa (x,y) = (d^d^5 aX + d { x x) d^) G (x,y) (5.11) 
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The Green's function of the half plane Go(x, y) is defined in terms of g-y and g 2 , the real parts 
of the analytic functions ln[(xi —yi) + i(x 2 —y2)] and ln[(— x\ —yi) + i(x 2 —1/2)] respectively. 
The imaginary parts of these functions are y) = -j^ arctan( ^~^ ) and $ 2 (x,y) = 

^ arctan( _f^ i ~^ i ) respectively. Therefore we can use the Cauchy-Riemann equations to get 

dWg 1 (x,y)=e,„dM$ 1 (x,y) (5.12) 

and 

dWg 2 (x,y) = -e»„dW<S> 2 (x,y). (5.13) 
With these relations we write the partition function as 

Z\Ap, J]=K. exp (-^j: J n d2x d 2 y A^t ^{x, y)A v (y) 



exp (^f Q d2 x d 2 V [B(x)$ 1 (x,y)J(y) + J(x)<f> 1 (x,y)B(y)} 
exp Q f Q d 2 x d 2 y [B(x)$ 2 (x,y)J(y) - J(x)^ 2 (x,y)B(y)}^ 
exp (| Jj 2 x d 2 y J(x)G (x,y)J(y)^ (5.14) 

We can use again the fact that B(x) = e^d^A^x) to obtain 

Z\Ap, J] = /Cexp {^J Q d 2 x d 2 y [^B(x)g 1 (x,y)B(y) + nJ(x)g 1 (x,y)J(y)] 
exp {^J Q d 2 x d 2 y [iB(x)$ 1 (x,y)J(y) + iJ(x)<S> 1 (x,y)B(y)} 
exp (\ f d 2 x d 2 y [-B(x)g 2 (x,y)B(y) - tt J (x) g 2 (x , y) J (y)\ 
exp Q J^d 2 x d 2 y [iB(x)$ 2 (x,y)J(y) - tJ(x)<S> 2 (x,y)B(y)}^ . (5.15) 

Replacing J{z) = a5(z — x) + b8(z — y) and B(z) = 7i(—5(z — x) + S(z — y)) in the 
integrands, the partition function 

Z\Ap, J] = /Cexp (7r(— 1 + ab)gi(x, y) + m{—b + a)$i(x, y) — 27r(& + a)$ 2 (a;, ?/)) 

(7T 7T \ 

-(1 - a 2 )^ 2 (a;, x) + -(1 - b 2 )g 2 (y, y) - tt(1 + a% 2 (z, y).J (5.16) 

Then, 
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Z\A , J] 

< exp (iau)(x) + ibuj(yj) >' ~~ 



feos Z[0] 

= exp(7r(-l +06)5-1 (a;, y) +iic(—b + a)$ 1 (x 1 y)) 
exp (— 7r(l + ab)g 2 (x, y) - in(b + a)$ 2 (x, y)) 



exp 



Q(l - a 2 )^ 2 (x,x) + ^(1 - & 2 )<? 2 (?/,?/)) (5.17) 



There are four posible combinations of a and b, which correspond to the four matrix elements 
of S (x,y). If a = 1, b = —1 we get 



(exp (icj(x) — iu {y))) = exp 



In ((x 2 - ?/ 2 ) 2 + (xi -2/1) 2 )' + 



+z arctan(— — — + m) 
xi - y% 



-. (5.18) 

z 



Comparing with Eq. Q5.4| ), we see that this term corresponds, in the fermionic theory, to 
Sq 2 (x, y) =< R(x)R'(y) >. In the same way, with a = b = — 1 we get 



(exp (—iuj(x) — iui(y))) = exp 



In ((x 2 - y 2 ) 2 + (-£! -yi) 2 )' 2 + 



+2 arctan 



1 

w 



(5.19) 



K -x 1 - 

which corresponds to S 22 (x,y) =< L(x)R'(y) >. 
Hence, we can identify 

R(x)~ exp(+iu(x)) JC- n (x) 
R*(x)~ exp(—iuj(x)) /C +7r (x) 
L(x)~ exp(— ioj(x)) /C +7r (x) 

L f (x)~ exp(+iw(x)) /C_^(x) (5.20) 

where ~ means that we identify products of operators (in the sense of the operator product 
expansion) inside arbitrary matrix elements of the fermionic and bosonic theories respec- 
tively. In the same way we can identify the remaining matrix elements. We conclude that, 
for the case of the boundary condition R = —L, the fermion operator on the half-plane is 
constructed in the same way as in FK [F2^]. The boundary condition appears only through 
the presence of diagonal operators which mix L and R components of the Fermi field. 
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It is worthwhile to note that, under global chiral (euclidean) transformations, R — > 
e e ° R and L — ► e -6 * L, the off-diagonal matrix elements of the fermion propagator remain 
invariant but the diagonal matrix elements do not. In the bosonized theory, global chiral 
transformations are represented by to — > uj + 6q. This representation of the matrix elements 
clearly satisfies this symmetry. Notice, however, that a global chiral transformation requires 
that the boundary condition be modified since the boundary condition itself breaks the 
chiral symmetry explicitly. Below we will come back on this issue. 



B. Dynamical Boundary Conditions 

Up to this point, we have shown how, following FK, we can calculate the one particle 
Green's function in its bosonized form for the boundary conditions R{x2) = —L(x2). The 
next step is the calculation of the Green's function with dynamical boundary conditions. In 
other words, we are looking for the solutions of 

x S F (x,y\e) = 5(x-y) 
S F \0,x 2 ;y\6) = -e 2e ^Sf(0,x 2 ;y\9). 
Here, 6{x 2 ) represents a quantum mechanical degree of freedom localized at the boundary, 
i. e. a quantum impurity. For the rest of this section we will take 9{x 2 ) to be arbitrary 
but fixed. In the next section we will consider a model with a fully dynamical, quantum 
mechanical boundary degree of freedom. 
The solution of Eq. ( |5.21| ) is 

/ b{x)a'\y) b{x)b'\y) \ 



S F (x,y\0) = — 



w 



a(x)a 1 (y) a(x)b x (y) 

y z w J 



(5.22) 



The functions a(x) and b(x) have to be chosen in such a way that the differential equation and 
the boundary condition of Eq. ( |5.21| ) are satisfied. In particular, it has to be (d X2 +id Xl )b(x) = 



(d X2 - id Xl )a(x) = 0. 

To find the explicit form of a(x) and b(x) we will use the FK approach. In the 
method developed in the previous subsection, the bosonic field satisfies vanishing bound- 
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ary conditions. This reflects the fact that the boundary conditions for the fermions were 
R(0,X2) = — L(0, x 2 )- But we have yet to determine how dynamical boundary conditions 
may affect the boundary conditions for the bosons. Our strategy will be to relate Sf(x, y\9) 
with a Green's function that satisfies non-twisted, R = —L, boundary conditions. For that 
purpose we define S(x, y) as a solution of the following equation 

(i$ x +fi)S{x,y)=5{x-y) 
S u (0,x 2 ;y) = -S 21 (0,x 2 ;y) 

Now, it is obvious that we can find S(x, y) by using the method described in the previous 
subsection. To establish a relation between S(x,y) and Sp(x,y\9) we write as 

4 = e ^P ( 5 - 24 ) 

where p(x) is a function such that p(0, x 2 ) = 9(x 2 ). It can be seen that if we take, 

S F ( X , y\6) = e t5p{x) S{x, y)e^ p{y) . (5.25) 

then Sp(x,y\9) is a solution of Eq. ( p.21| ). Following the same steps that lead to the calcu- 
lation of So, we define 

Z[Afj,, J, Afj\ = J Vuj exp S b (uj, A^) 

exp (i jf J(z)cu(z) + J Jp(z)Ap(z)\ 

exp (— J dx 2 dy 2 9(x 2 )K(x 2 ,y 2 )9(y 2 ) + 

^ J dx 2 dy 2 A 2 (x 2 )K(x 2 ,y 2 )9{y 2 ) s j (5.26) 

Here, A^ and J are the external sources defined in Eqs. ( |5.6| ) and Q5.8| ) respectively, and 
Jfi(z) is the current defined in (|4.9|) . If we define Sb(9, A 2 ) as 

S B (9,A 2 ) = ^- J dx 2 dy 2 9(x 2 )K(x 2 ,y 2 )9(y 2 ) + - J dx 2 dy 2 A 2 (x 2 )K(x 2 ,y 2 )9(y 2 ) (5.27) 

the partition function can be written as 
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Z[A„ J, A M ]= e s »^ J Vu exp (s B (u, + i JjJ(z) - ^d u A^z)]cu(z) 

= e SBie ' Aa) J Vu exp (s B + iJ J(z)u(zj\ (5.28) 

with J{z) = J(z) — -e^dyA^z). This expression is formally equal to Eq. (p.7|). Hence, we 
can follow the same steps that took this equation into Eq. (|5.10|) to obtain 

Z[A^, J, A^]= e 3 -^ Z[A„ J] 

exp Jjd,A,(x)G (x,y)V 2 p(y) + V 2 p(x)Go(x,y)d,A,(y)]) 

exp |jJ(x)G (x,y)V 2 p(y) + V 2 p{x)G {x,y)J(y)] S j 

exp (i- 1 4(x)r^(x, y)A v {y)j (5.29) 

where Z\A^ J] is exactly equal to Eq. ( |5.15| ) and the other factors come from the extra term 
in J(x). These two factors can be computed with the same techniques used before giving 

exp (— ap(x) — bp{y)) exp ( — / dz 2 p(Q,z 2 ' ' ' " ~~ M 



+ (x 2 - z 2 ) 2 J 

exp ( - / dz 2 p(0, 22) T , / V * Z — ^ I • (5.30) 

2/i + (2/2 - z 2 y J 



s 7T 

7T 



Then, 



(exp (iau;(x) + ibu(y)))i 



Z[A^ J, A^j 



2^] 

= exp (ir(— 1 + ab)gi{x, y) + in(—b + a)$i(x, y)) 

exp (— 7r(l + ab)g 2 (x, y) - 27r(5 + a)$ 2 (^, J/)) 
1 /" , / n * a^i + i(x 2 - z 2 )\ 



/ 1 r axi + i[x 2 - z 
exp - / cte 2 p(0, 22) ; 

\7l J Xf + [X 2 - Z 2/ 

f 1 /" , , n N 6yi - %2 - 22) 
exp I — y cte 2 p(0, z 2 



J/l + (2/2 - Z 2 f ) 

exp (— ap(x) — bp(y)) (5.31) 

In this case also, there are four combinations of a and b which give the four matrix elements 
of S(x,y). Up to a normalization constant the one particle Green's function is 
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where 



S(x,y) = - 



( e-P^b(x)g- 1 (y)e-P^ e -P^b(x)b' 1 {y)eP^ \ 
w z 

eP ( - x )g{x)a- 1 {y)e~PM e P( x ~> ajx^- 1 (y)ePW 
z w 



(5.32) 



(5.33) 



and 



6 W =exp(i/^9(z 2 ) V y l + + '^_ z f 



Note that p(x\,x 2 ) satisfies 



It can be seen that 



V 2 p{x) = in xi > 
p(0, x 2 ) = 9{x 2 ) 

\na(x) = -2[ ^ p (^)e^ 2 -^ 
Jui<o 2n 

= -2p_(rr 2 - ixi) 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



Hence, lna(x) is the anti-analytic extension of the negative frequency part of p(0,x 2 ) 
9(x 2 ). In the same way, 



\nb(x) = 2[ ^p{u)t 

Jw>0 2,71 



,iu}(x2+ixi) 



2p + (x 2 + lXi) 



(5.38) 



is the analytic extension of the positive frequency part of p(0, x 2 ) = 6{x 2 ). Note that 
p(x\,x 2 ) is the real part of both, p+(xi + ix 2 ) and p~(xi — ix 2 ) and that 



2p + (x 2 + ixi) + 2p_(rr 2 - ixi) = 2p(xi ) x 2 ) 



(5.39) 
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Using the relation between Sp(x, y\9) and S(x,y) (Eq. ( |5.25|) ), we get an expression in 
the bosonized theory for the fermion Green's function which satisfies dynamical boundary 
conditions (Eq. ( |5.22| )). 

By using Eqs.( |5.33| - |5.34| ) we can rewrite Eq.( |5.22| ) as 



S F (x,y\6) 



2tt 



( e *J dz 2 0(z 2 )h 11 (z 2 ) g l J dz 2 8(z 2 )h 12 (z 2 ) ^ 

W z 
A j d.z 2 0(z 2 )h 21 (z 2 ) i j dz 2 g(z 2 )h 22 (z 2 ) 



(5.40) 



where 



foil (22) 
hi 2 (z 2 ) 
h 2 i(z 2 ) 

h 22 (z 2 ) 



1 



+ 



1 



xi - i(x 2 - z 2 ) y% + i(y 2 - z 2 ) 
1 1 



xi - i(x 2 - z 2 ) y 1 - i(y 2 - z 2 

1 1 

+ 



xi+^xa-^) yi + i(y2-z 2 ) 

1 1 



(5.41) 



xx+ifa-z^ yi-i(y 2 -z 2 ) 

We can summarize these results in the form of a set of bosonization rules for the fermion 
operators, which now include the effects of the presence of a boundary degree of freedom. 
As in the previous case, i. e. for R(x 2 ) = —L(x 2 ), we can identify the fermionic operators, 
R(x 2 ) and L(x 2 ), with the corresponding bosonic operators. From Eq. (|5.22 ) we see that 

1 



-R(x 2 )~ exp iu(x) + 
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R\x 2 )~ 
L(x 2 )~ 



1 



exp 
exp 



-IUJ{X) 

7T 



-iu)(x) 



71 



L(x 2 )~ exp +iu(x) 

\ ' 7T 



dz 2 9 
dz 2 
dz 2 
dz 2 



, , Xi + 1{X 2 - z 2 
Z *)—2 



*2. 



^2, 



xf + \x 2 - z 2 y j 

x\ + i(x 2 - z 2 ) 



x\ + (x 2 - 




—x 1 + i(x 2 




%i + (x 2 - 


z 2 y 


—xi + i(x 2 


-Z2) 



x\ + (x 2 - z 2 ) 



IC +n (x) 
/C +7r (.x) 
IC-Jx) 



(5.42) 



where, as before, ~ means that the bosonic operators yield the same correlation functions 
(with the same boundary conditions) as the operators R and L in the Fermi theory. This 
result completes the bosonization construction. Notice that the main modification is the 
presence of a boundary operator m the definition of the fermion operator of this generalization 
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of Mandelstam's construction. We should stress here that, as remarked in |T|, in order to 
compare with the operator construction in Minkowski space, not only the imaginary time 
x 2 needs to be Wick-rotated back to the real time x$ but the boundary degree of freedom 9 
also has to be continued to the imaginary axis, 9 — > i9. 

VI. FREE FERMIONS COUPLED TO A DYNAMICAL BOUNDARY DEGREE 

OF FREEDOM 

Now we will compute the one particle Green's function in the case that 9 is a quantum 
mechanical degree of freedom. That is, instead of regarding 9(x 2 ) as a prescribed classical 
parameter, it has to be regarded of as an extra degree of freedom of the full system. In order 
to account for the effects of the dynamics of 9, we have to specify its dynamics through an 
extra term in the lagrangian and we must also integrate over all possible configurations of 
9. Hence, the full one particle Green's function becomes 

JVi>VipV9 tp{x)ip(y) e iS$,M)+ism 



Sf< ^ x ^^ j VipV%l)V9e iS ^> e )+ iS ° (») 

JV9 S F (x,y\9)Det(9)e iS ^ 

(6.1) 



JV9 Det(0) e iS ^ 

where Sty, ip, 9) is the action for the masless fermions in the half line coupled to a dynamical 
degree of freedom at the boundary. Det(#) is the determinant calculated in Section III 
(Eq. ( |3.39| )). Upon exponentiation, it yields an induced term in the action S in d(9) of the 
form 

( V>et{i$) m -i \ 1 r [9(y 2 ) - 9{x 2 )f 

S ind {9) = In -— w = — / dx 2 dy 2 - — 6.2 

\ Uet(^) B J 4:ic 2 J {x 2 -y 2 ) 2 + a 2 

which represents the dynamics of 9 induced by the fermionic degree of freedom. The intrinsic 
dynamics of the boundary degree of freedom is represented by the term Se(9) in the action. 

At this point, it is convenient (and helpful to understand the physics) to make an analytic 
continuation of the boundary degree of freedom 9 — > i9. In this way, we recover the usual 
periodicity property 9 — > 9 + 2tt . In Euclidean space this symmetry is not manifest since 
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the group of chiral transformations is non-compact in Euclidean space. This change has 
for effect to (a) change the sign of the induced action S in d{9) and (b) the Green function 
amplitudes need to be analytically continued (Eq.( |5.40| )). 

The intrinsic dynamics of the boundary degree of freedom 9 is determined by the un- 
derlying microscopic system from which the Luttinger-type model is derived. A simple and 
still realistic model of this dynamics consists in regarding 9 as a coordinate for a single 



S e {9)= jdx 2 \^{^-\ -V{9)\ (6.3) 



degree of freedom moving on a ring. This is the situation in the Caldeira-Leggett model |30 
which describes the effects of a local degree of feedom coupled to a "fermion bath". From 
the point of view of a Hubbard model or a quantum wire, 9 represents a boundary effect 
of electron-electron correlations. A set of "interesting" terms that can be incorporated in 
Sg(9) are (in imaginary time) 

, V 

2 \dx 2 

where M is the mass of the boundary degree of freedom. In quantum wire problems, the 
potential V (9) is the backscattering amplitude at the boundary. This term is of the form 

V(6) = Gcoa(J36) (6.4) 

where (3 is a parameter which, if interactions are present, depends on the electron-electron 
coupling constant. 

By simple inspection we can see that all terms of the action, except for the potential 
term, are invariant under a global chiral rotation of the boundary conditions, 9(x 2 ) — > 
9{x 2 ) + 9 Q . The potential term reduces this continuous symmetry to a discrete subgroup 
9(x 2 ) — > 9(x 2 ) + 2irn//3 , where n is an arbitrary integer. Similarly, the off- diagonal terms 
of the Green function are invariant under global chiral rotations whereas the diagonal terms 
are not. These symmetries have a very simple interpretation in the bosonized theory. Since 
the chiral symmetry of the fermions translates into the invariance of the bosonized theory 
by constant shifts of the field u, the boundary condition u = at X\ = reflects the fermion 
boundary condition R = —L. The boundary degree of freedom 9 can be viewed as the 
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value of uj at the boundary. Hence, a path integral in which one integrates over all values 
of 9 without restrictions should correspond to free boundary conditions for u. We will see 
below that this is indeed the case and that, in the absence of potential terms, the boundary 
conditions are free and the system behaves as in the case of an infinite line. Conversely, 
potential terms act to "pin" the boundary value of the bosonic field u to a finite set of 
possible values. This behavior is closely related to the picture of Affleck and Ludwig of the 
Kondo problem as a crossover in boundary conditions [ 17| . 



In the framework of macroscopic quantum tunneling and coherence, Fisher and Zw- 



erger |3T| have studied systems with actions of this same type. By means of renormalization 
group arguments, they concluded that the kinetic energy term is irrelevant at low frequen- 
cies (or long times), that the induced term in the action is a strictly marginal operator 
and that the cosine terms of the potential are relevant. Fisher and Zwerger refer to this 
strong-coupling fixed point as to the localization transition. The induced term in the action 
plays the same role in this problem as the free massless scalar theory (or Gaussian model) 
does in two-dimensional critical systems and in conformal field theory. 

We are interested in the computation of the matrix elements of the fermion one-particle 
Green's function for a system with dynamical boundary conditions. To simplify the discus- 
sion we will include the effects of the potential term only within a linearized approximation 
in which V(9) ~ — ^-9 2 = — ^9 2 . This approximation misses the important periodicity 
property and overemphasizes the physics of localization. Recent work by Ludwig and colab- 
orators and by Tsvelik indicates that there is interesting physics which this fixed 
point does not describe correctly. 

Hence, we conclude that the effective dynamics of the boundary degree of freedom 9 is 



described by the action 

1 c c ( A/f / f)9 \ 2 C \ 

S eS (9) = —j dx 2 dy 2 9(x 2 )K(x 2 ,y 2 )9(y 2 ) + J dx 2 I — f — J + -9 2 \ (6.5) 

where K(x 2 ,y 2 ) is given by Eq.( |4.8| ). In order to compute the one-particle fermion Green's 
function we need to evaluate functional integrals over 9 of the form 
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W a , = (exp ( - J dz 2 e{z 2 )K,{z 2 ))) 6 = JV6 e W (-S eS (6)) (6 ' 6) 



Since the exponent is quadratic in 9, these integrals can be done explicitly. In Fourier space, 
the action S e s(9) is 

SA») = J^)\^ (6.7) 
where the kernel K e ^{uj) is given by 

K eS (u) = l ~ e ' ' + ttMuj 2 + 7iG (6.8) 
a 

The form of the kernel K c g(u) shows that the inertial term acts like a high frequency cutoff. 
At low frequencies we can use the approximation K eS (uj) ~ |u;| + nG. The expectation 
values of the form of Eq. ( |6lf ) become 

1 r cL< , 1 

^eff 



W af3 = exp(--^- f ^h a p{-u) 1 , Kp{u)) (6.9) 

Z7T J ZlX Kpff(Uj) 



where h a p(uj) is the Fourier transform of h a p(z 2 ). By making use of the expressions for the 
amplitudes h a p{z 2 ) given in Eq.( |5.41|) , we find that their Fourier transforms can be written 
in terms of the functions 

/+oo pivr 
dr - = 2ne-^ [uTiv)) ®{±uj) (6.10) 

-oo u±i(r — v) 

where u > and 0(x) is the step function. The Fourier transforms become 

h xl {uj) = F + {u',x u x 2 ) + F-(u;y 1 ,yz) = 2ne-^ {xi - iX2) Q(u) + 27ce^ liyi+iy2) e(-uj) 

h 12 (u) = F+(u;x u x 2 ) - F + (u J ;y 1 ,y 2 ) = 2vr ( e -M - e -M(w-iw)) e(uj) 

h 21 (uj) = -F-(u;x 1 ,x 2 ) + F-(cj;y 1 ,y 2 ) = 2tt (- e -H(*i+*«) + e -M(w-Nw)) e(-u) 

h 22 {u) = -F_(uj;x 1} x 2 )-F + (uj;y ll y 2 ) = -27re-^ +iX2 ^Q(-u) - 2ne-^- iy2 ^Q(uj) (6.11) 

We can now write down an explicit formula for the logarithm of the amplitudes W a p. We 
find that the off-diagonal matrix elements vanish identically, In W\ 2 = In W 2 \ = and that 
the diagonal matrix elements are given by 
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oo e -w((a!i+yi)-i(xa-jte)) 



\nW = \nW u = \nW 22 = -2 did — (6.12) 

Jo iv e ff(o;) 



The explicit computation of the integral in eq. 6.12 yields the result 



In W = 2e nG ™Ei(-nGw) (6. 13) 

where Ei(—x) is the exponential-integral function. This result yields the asymptotic behav- 
iors of In W as 

lnW«2(C + ln(-w7rG)), for \w\ttG < 1 

\nW & |— , for|w|7rG>l (6.14) 

nGw 

where C is Euler's constant. In the limit x\ + y\ — > 0, these results hold provided that 



X\ + yi is replaced by the smallest of the cutoff a and V ttM. 

With these results at hand we find that, for dynamical boundary conditions, the matrix 
elements of the Fermion Green's function are given by 

/ w 1 \ 



Sf{x,v) = 

Z7T 



z 



1 W 
\ z id J 



(6.15) 



where W is given in Eq.( |6.12|) . 

Hence, we find that the off diagonal matrix elements of the fermion Green's function, 
i. e. the propagators (R(x)Iv(y)) and (L(x)L'(y)) respectively, are equal to their value in 
the infinite plane and are unaffected by dynamical boundary conditions. Notice, however, 
that for a specified configuration of the boundary chiral angle #(#2), these matrix elements 
are modified. In contrast, we find that the diagonal matrix elements, i. e. the propagators 
(R(x)L^(y)) and (L(x)R Jf (y)) respectively, are modified by dynamical boundary conditions. 
However, these changes are only significant either close to the boundary and at short times. 
At long times and far away from the boundary, the off diagonal matrix elements behave 
exactly as in the case of R = —L boundary conditions. Close to the boundary, we find 
power law corrections of the form 
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w+ . const. (irGw) 2 
(R{x)L\y)) A- — i- (6.16) 

If the potential term is absent (G — > 0), the boundary degree of freedom 9 is unpinned. 
In this regime, the diagonal matrix elements of the fermion one-particle Green's function 
vanish identically, even away from the boundary, 



S F {x,y) = — 



v i / 

\ 2 / 



(6.17) 



In this limit the boundary degree of freedom fluctuates wildly and it effectively removes the 
boundary condition uj = for the bosonized theory. In other terms, the bosonized theory has 
free boundary conditions instead of Dirichlet. In this regime, the fermions are not scattered 
by the boundary. In a sense, they are "eaten" by the boundary degree of freedom. The 
fermion Green's function in this limit is equal to the Green's function in the infinite plane 
(1 e. the full line). Once again, this result agrees with the picture of Affleck and Ludwig. 

Finally, let us note that the amplitude W has a clear physical meaning. If W was 
a constant, it would have the straightforward physical interpretation as a backs cattering 
amplitude. However, it is both position and time- dependent. Hence, the backscattering 
amplitude determined from W is both momentum and energy (frequency) dependent. The 
actual form of this dependence is rather involved, as can be seen by inspecting the form of 
W in real time. The frequency dependence of W is simply a manifestation of the crossover 
between the two asymptotic behaviors shown in eq. |6.14| . Indeed, at frequencies high com- 
pared with 1/(ttG), W scales to zero at low frequancies, i. e. it behaves asymptotically as in 
the case in which the boundary degree of freedom is unpinned (G) — > 0). In this regime the 
effective value of the diagonal matrix elements of the fermion Green's function scale to zero 
and the effective boundary condition is free. Conversely, at low frequencies and for G fixed, 
W scales to a constant. This is the behavior of the pinned boundary degree of freedom and 
it amounts to a boundary condition R = —L for the fermions. 
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VII. CONCLUSIONS 



In this paper we have reexamined the bosonization of a theory of free fermions in the 
presence of boundaries and with specific boundary conditions. This is a very simple system in 
which the role of boundary conditions and of boundary degrees of freedom can be investigated 
quite explicitly 

We have exhibited the role of the boundary conditions of the fermions in the bosonized 
theory. We used the Forman's method for the computation of fermion determinants to 
derive the action of the bosonized theory which exhibits the effects of the boundary degrees 
of freedom. This action allowed us to determine the form of the currents and found that 
the currents acquire additive corrections due to the boundary degree of freedom. These 
corrections are needed in order to to satisfy the conservation laws. In contrast, we find that 
the current algebra is not affected by the presence of boundaries and of boundary degrees of 
freedom. We also constructed the bosonized form of the Fermi operators . We showed that 
the boundary degrees of freedom enter explicitly in the definition of the Fermi operators. 
Finally, we used these results to calculate the fermion one-particle Green's function for a 
theory with a dynamical boundary degree of freedom. 

The methods described in this paper can be generalized so as to include the effects 
of interactions and of non-abelian symmetries. It is also possible to use these methods 
to investigate the interesting non-perturbative structure recently found , even in abelian 
systems, by means of the thermodynamic Bethe Ansatz p2| , |33l . In a separate publication [fjj 



we apply our methods to the case of interacting fermions, in particular to the Luttinger model 
and study the interplay between the Kondo effect and interactions in 1 + 1-dimensional 
systems. 
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APPENDIX A: BOUNDARY JACOBIAN 



The aim of this appendix is the computation of 

\nSJ F = lim ( [ dt [ dx 2 [ dx\ < x\e~ s ^ >t ^ >t 

s,e^0 \Jo J Jo 



X > 4>{Xl, X2 



-2 lim 



dt / dx 



■ 2 j dxitr[K 2 (ip t ) x, x)j 5 (f>(xi, x 2 )} ) . 



(Al) 



In order to do this computation we have to find the heat kernel K 2 in (0, e). The operator 
ip can be written as 



ip 



-dt+B 
d x + B 



(A2) 



with B = id 2 + A 2 . Following APS [2l| we make the aproximation A 2 (x\, x 2 )~A(0, x 2 ). In 
order to compute the trace, we need to find a basis of eingestates of ipt- We construct such 
a basis as follows. Let us call e\(x 2 ) the eigenvector of eigenvalue A of the operator B. That 



is. 



Be x (x 2 ) = \e\(x 2 ) 



(A3) 



It can be shown that there exist w k ,x and a set of functions f k ,x(xi) and g k ,x(xi) such that 



ip^ K x{xi,x 2 ) = w k ,xipk,x(xi,x 2 ) 
fk,x(0) = -9k,\(P) 
where by definition ipk,x( x ii x 2 ) is 



^k,\{xi,x 2 ) 



( fk,x(xi)e x (x 2 ) ^ 



gk,x{xi)e x {x 2 ) 

^ R k ,x(xi,x 2 ) 
L ktX (xi,x 2 ) 



J 
\ 

) 



Specifically, 

^k,x{xi,x 2 ) 



e x {x 2 ) 



(2(A + w k , x )w kt x) 2 



( 



(A + Wk t x) sin kx\ — k cos kxi 
(A + Wf.,x) sin kxi + k cos kx\ 



\ 
) 



(A4) 



(A5) 



(A6) 
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form the complete set that satisfies Eq. ( |A4j ) if we take = (k 2 + A 2 ) 3 (k > 0), and 
A = Y~(n + | + A) where A is the flux of the gauge field through the boundary. 
Then, 



\n5J F = lim (-2 I dx 2 (j)(0,x 2 ) f dx x [Y r! AxMe' 8 ^^ R k x (x 

£4>i)0e-^+^,A(*i)] 

k,X j 

r r°° dk ( P - s i k2 + x2 ) \ 

= lim / <fe 2 0(O, * 2 ) / - £ — — r cos(2fce - 1). (A7) 

It is important to note that in the above limit, s < e. In the limit T —>■ oo the sum over A 

can be replaced by an integral. By using the Euler's formula it can be shown that 

—2T r r°° dk r°° ( e - s ( k2 + u2 ) \ 
\nSJ F = — — lim / dx 2 0(0, x 2 ) / — / du\ r cos(2A;e - 1) (A8) 

TY 2 s,e^0j J 7T Jo \(k 2 +U 2 )2/ 



The result of this last integral is 



-^i 1 -*^)} <A9) 
where 5 is a Schwinger parameter that goes to zero and is the error function. The 

order of the limits is 5 <C s <C e. Then in the first term of Eq. ( |A5| ) we use the asymptotic 

expansion of $(a/x). But in the second term we use the limit for z C of Hence 

Eq. ( |A9| ) becomes 

\nSJ F = — lim / dx 2 0(0, a; 2 )^ f F ^ )v ^ - . (A10) 



47T s,e-*0i ' y^S 

The second term in Eq. ( |A10|) tells us that there is a singularity in s. According to the 
definition of the heat kernel regularization p] we have to keep only the finite part of ( A10 ). 
Therefore, the properly regularized jacobian at the boundary is 

SJ F = exp (-— J dx 2 0(0, x 2 ) \ . (All) 

Since 0(0, x 2 ) = 0, the contribution of this jacobian to the partition function is one. 
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APPENDIX B: FORMAN'S DETERMINANT 



In this appendix we are going to set up the notation and and show some of the calculations 
that lead to eq. ( |3.35| ) First note that we are working in the space spanned by e lWnX ° with n 
runing over the integers. This space splits naturally into the one expanded by e lWnX ° with 
n > (the space of functions with positive frequencies) and the space spanned by e lWnX ° 
with n < 0. Indeed, each a can be decomposed into a + and ct_ such that a = a + + a_ and 



n>0 

a_(:r 2 ) = 5> n e^ 2 . 

n<0 



(Bl) 



In the same spirit, we can decompose operators A acting in such space into blocks: 



A 



( \A\__ [A)_ + ^ 



[B2) 



y[A) + _ [A] ++ j 

where [A] means [-A] m _ n with m, n < and so on. 

Finally, we use the notation a = a(/i, x 2 ) = fia(x 2 ) The first step is the calculation of 
the inverse of 

\ 



h(n) 



( [e- c 



\ 



(B3) 



Formally it has the form 



( [U}__ [U]_ + ^ 



(B4) 



where 



,-i 



[U\- = {1 " ([e- ot ]„)" 1 [e a ]_ + ([e fl '] ++ )- 1 [e-«] + _}" 1 ([e-»]_)- 
[U} ++ = {1 " ([e a ] ++ )- 1 [0 + _([e- a ]_J" 1 [e1_ + } -1 ([e a ] ++ )- 1 , 

PU = ([c- a ]__)" 1 [c a ]^{([c a ] 44 .)- 1 [c- a ] + _([c- a ]__)~ 1 [c a ]_ + " 

[U]+- = ([e a ] ++ )~ 1 [e- a ] + _{([e- a ]__)- 1 [e Q ]_ + ([e a ] ++ )- 1 [e- a ] + _ - 1}'%-°]. 
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(B5) 



We need to compute 

{1 - ([e-°]__)" 1 [e"]_ + ([e°] ++ )- 1 ([e-°] + _)}" 1 ([e-'-]__)- 1 } 

+ I >{(|;« e i ++ )->-i + -)«^-r I M- + ) 

{1 - ([e1 ++ )- 1 [e-1 + _([e-«]__)- 1 ([e«]_ + )}" 1 ([e«] ++ )- 1 } . (B6) 

In order to compute this last expression it is usefull to derive some identities relating the 
matrix elements used above. Note that for deriving this expressions we are going to use the 
fact that a is an abelian variable. Appart from that a could be any antiperiodic function in 
[0, 0\. The main propertie that we are using here is the fact that [e a ~] p = for p > and 
[e a +] n = for n < 0. Having this in mind the following identities are straightforward 





l + - = o 


(B7) 


[ e «+; 


U = o 


(B8) 


[e a -\ ++ [e a+ ; 




(B9) 


[e a -\ ++ [e- a -\ 




(BIO) 


[e a+ } ++ [e- a+ [ 


! ++ = i++ 


(Bll) 


[e«+]__[e«-; 


|__ = [ e a +e Q -]__ 


(B12) 




1- = 1- 


(B13) 


[ e °-]__[ e - a -; 


1- = 1- 


(B14) 




l ++ = [e a -e° + ]_ + 


(B15) 


[e^]__[e a -; 


U = [e a+ e a -U 


(B16) 


[e a -] ++ [e a \ 


\ + _ = [e*-e Q +] + _ 


(B17) 


[ e ^] + _[ e °-; 


|__ = [e a+ e Q -] + _ 


(B18) 
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[e a+ ] ++ [e~ Q+ ] + _ = -[e a+ ] + -[e- a + ]__ (B19) 
[e a -]__[e- a -]_ + = -[e a -}_ + [e-«-} ++ . (B20) 

With a change a — * —a more identities can be derived. It is inmediate to realize that 

([e-U)- 1 = [e«-]__[e«+]__ 
([e1 ++ )^ = [e^ + ] ++ [e- Q -] ++ . (B21) 

As an example of how to use the identities (B7)-(B20) we compute the following expression 
{1 - ([c-]„)- 1 M_ + ([c-] ++ )- 1 ([c-] + _)}- 1 . 
We define [M\ as 

[M]__ ee ([ e - a ]__)- 1 [ e a ]_ + ([e a ] ++ )- 1 [ e - a ] + _. (B22) 
By using Eq.(B21) we obtain 

[M\__ = [e a 1„_[e« + ]„[e1_ + [ e -« + ] ++ [ e — ] ++ [e-] + _ (B23) 

With Eqs.(B15) and (B18) we can write the factors of the form [e a ] in terms of [e a ~] and 
[e«+] 

[M]__ = [e a -]__[e Q +]__[e Q -]_ + [e^] ++ [e- Q + ] ++ [e- Q -] ++ [e- Q + ] + _[e- Q -]__ = 
= [e a -]__[e a +]__[e a -]_ + [e- a - ] ++ [e-+] + Je— ]__. 

After rearranging some factors using Eqs.(B16) and (B17) [M] can be written as 

follows 

[M]__ = [e«-]_Je a ]_ + {e- a ] + Je- a -]__. 

Hence 

(1 - [M]_J = [ e °-]__(l - [e Q ]_ + [e- Q ] + _)[e- a -]__ 

= [e a -]__[e Q ]__[e- a ]__[e- a -]__. (B24) 

By using eqs.(B14) the inverse of the above expression reads 
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(1 - [M]__y 1 = [e— ]__([e-1_J- 1 ([e«]_J- 1 [e-]__. (B25) 

The factors ([e~ a ]__y 1 and ([e a ]__) _1 can be computed with the help of eqs.(B12)-(B14). 
The result is 

([e- a ]__)- 1 ([e a ]__)- 1 = [e-]__[e^]__[e— }__[e—}__. (B26) 

Therefore, the expression we were calculating turns out to be 

{l-([e- Q ]__)- 1 [e«]_ + ([e Q ] ++ )- 1 [e-«] + _} = 

= [e Q -]__[e a -]__[e a +]__[e- a -]__[e- Q +]__[e- Q -]__. (B27) 

Following the same steps we can see that 

{1 - ([e ffl ] + J- V^IO- J~V]- + } = 

= [e- a+ ] ++ [e- Q+ ] ++ [e- a 1 ++ [e a+ ] ++ [e a -] ++ [e a+ ] ++ . (B28) 

The next step is the calculation of the derivatives. 

(A[ e «] n _J = -J e ~ 2ni{n ~ m)x2 a(x 2 )e- fia{x2) 



= -E 



-/ia(x )l 



ln-m-pi a \p 



-J2[e-^] n _ q [a} q _ m . (B29) 
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And in our notation, 

'dp 



(4r.[^-) = -([e- Q ]__N- + [e- a ]_ + N + _). (B30) 



In the same way, 

d 



^ L e- a ]_ + = ([e1__N_ + + [e1_ + H ++ ) 
^[e1 + _ = -([e-1 + _N__ + [e- a ] ++ H + _) 



= ([e Q ] + _N_ + + [e°] ++ [a] ++ ). (B31) 



By inserting (B27)-(B31) into (B6) we obtain 
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Tr[j-h^)h~\^] = 

Tr [a]- + {[e- Q +] ++ [e-^] ++ [e—] ++ [e^] ++ [e«] + _ - 

[e- a+ ] ++ [e- a -] ++ [e-^] + _[e-]__[e^]__[e— }__[e a }__} + 
Tr [a] ++ {[e-^] ++ [e-^] ++ [e— ] ++ [e Q +] ++ [e*] ++ - 

[e-^] ++ [e— ] ++ [e-^] + _[e a -]__[e a +]__[e- a -]__[e a ]_ + } + 
Tr N + -{[e a -]_J e ^]__[e a -]_ + [ e -^] ++ [ e — ] ++ [e^] ++ [e- a ] ++ - 

[e°- ] [e a - ] __ [e a + ] [e— ] [ e -] + } 
Tr N-{[e a -]_J e ^]_J e a -]_ + [ e ^] ++ [ e —] ++ [ e ^] ++ [e- a ] + _ - 

[e a 1-[e a+ ]__n_Je- Q 1 ++ [e-l_}. 



(B32) 



By using the identities (B7)-(B20) we can see that the factor of [a] contributes with —1 

and the factor of [a\++ contibutes with 1 . The other two terms are more complicated to 
calculate but it can be done. The final result is 

Trl-^Hfijh- 1 ^)] = tr{[a] ++ - [a] — 

+ W + -(-[e Q -]__[e- Q -]_ + + [e a -]__[e a -]_ + [e- a -] ++ [e- a -] ++ ) 

+ H- + ([e- Q+ ] ++ [e a+ ] + _ - [e' a+ ] ++ [e- Q +] + _[e a +]__[e a +]__)} . (B33) 



APPENDIX C: EUCLIDEAN-MINKOWSKI CORRESPONDENCES FOR DIRAC 

FERMIONS 

Upon the analytic continuation from real time x$ to imaginary time X2 = ixo, a number 
of quantities are transformed. Firstly, the Minkowski space metric g^ v 



(CI) 



becomes, in euclidean space 
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(C2) 



The Minkowski space 7-matrices are (in the chiral basis) are 



70 = 0"! 7l = -«02 75 = 0"3 



(C3) 



and have the property 



tr(7 M 7^) = 2 g^u- 



(C4) 



The euclidean 7-matrices are chosen to be hermitean, in which case they are 



7l = °2 72 = Cl 75 = 0"3 



(C5) 



and they satisfy the relations 



tr(7 M 7„) = 2 V, 757^ = -ie^ v . 



(C6) 



We also give the correspondence for 



■0 = *0' 7o — > i*0 



(C7) 



as well as for vector fields, which transform like 



A -> iA 2 Ai -> Ai 



(C8) 



With this notation the covariant derivatives become 



d- z - iA z 



d z + iA s 



(C9) 



where we have set 



(CIO) 



d z — di - id 2 , dz — di + id 2 
A Z = A 1 - iA 2 , A- z = A 1 + iA 2 . 

Using these identifications the Minkowski space actions Sm and Se are related in the usual 
manner 
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iS M -> -S E . (Cll) 

For the specific case of a massive Dirac fermion coupled to a background gauge field, the 
Lagrangian Cm, 

L M = i> (i0+A)i^- Mvpip - iM^ls^ (Cf 2) 

becomes, in euclidean space, 

C E = i> (0 + A) + iM^ - M 5 ^7 5 ^. (Cf 3) 

Notice that the factor of % in the mass term results from our choice of hermitean euclidean 
7-matrices. Had we chosen antihermitean 7-matrices the factor of % would have been absent 
from the fermion mass term. 

Finally we note that the Minkowski space fermion propagator 

S$ = -i{f^ a {x) Mx')) (C14) 

(where (A) is the ground state expectation value of A and T is the time ordering operator) 
becomes, in euclidean space, 

Sf = M*'))- (C15) 
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